Differential Equations 



The subject of ordinary differential equations encompasses such a large field that you can make a 
profession of it. There are however a small number of techniques in the subject that you have to know. 
These are the ones that come up so often in physical systems that you need both the skills to use them 
and the intuition about what they will do. That small group of methods is what I'll concentrate on in 
this chapter. 



4.1 Linear Constant-Coefficient 

A differential equation such as 



d 2 x V 

w) 



+ t V + 1 = 



relating acceleration to position and time, is not one that I'm especially eager to solve, and one of the 
things that makes it difficult is that it is non-linear. This means that starting with two solutions Xi(t) 
and X2(t), the sum X\ + X2 is not a solution; look at all the cross-terms you get if you try to plug the 
sum into the equation and have to cube the sum of the second derivatives. Also if you multiply Xi(t) 
itself by 2 you no longer have a solution. 
An equation such as 

+ cfix . q dx 

may be a mess to solve, but if you have two solutions, X\(t) and Xi(t) then the sum ax\ + (3x2 is also 
a solution. Proof? Plug in: 

e , #(a^ + te) + ^(ax.+fe) _ {axi + 

This is called a linear, homogeneous equation because of this property. A similar-looking equation, 

^ d^x 2 dx 

e W + t dt~ x = t 

does not have this property, though it's close. It is called a linear, inhomogeneous equation. If X\(t) 
and X2(t) are solutions to this, then if I try their sum as a solution I get It = t, and that's no solution, 
but it misses working only because of the single term on the right, and that will make it not too far 
removed from the preceding case. 

One of the most common sorts of differential equations that you see is an especially simple one 
to solve. That's part of the reason it's so common. This is the linear, constant-coefficient, differential 
equation. If you have a mass tied to the end of a spring and the other end of the spring is fixed, the 
force applied to the mass by the spring is to a good approximation proportional to the distance that 
the mass has moved from its equilibrium position. 

— * 

If the coordinate x is measured from the mass's equilibrium position, the equation F = ma says 

hM«^D m § = -fa (4.1) 
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If there's friction (and there's always friction), the force has another term. Now how do you describe 
friction mathematically? The common model for dry friction is that the magnitude of the force is 
independent of the magnitude of the mass's velocity and opposite to the direction of the velocity. If 
you try to write that down in a compact mathematical form you get something like 

-* v 

-^friction = — /^k-Fkl -pzr (4.2) 

This is hard to work with. It can be done, but I'm going to do something different. (See problem 4.31 
however.) Wet friction is easier to handle mathematically because when you lubricate a surface, the 
friction becomes velocity dependent in a way that is, for low speeds, proportional to the velocity. 

^friction = ~bv (4.3) 

Neither of these two representations is a completely accurate description of the way friction works. 
That's far more complex than either of these simple models, but these approximations are good enough 
for many purposes and I'll settle for them. 

Assume "wet friction" and the differential equation for the motion of m is 

d x -, i dx , , 

m-r^ = -kx - b-r- (4.4) 
dt 2 dt 

This is a second order, linear, homogeneous differential equation, which simply means that the highest 
derivative present is the second, the sum of two solutions is a solution, and a constant multiple of a 
solution is a solution. That the coefficients are constants makes this an easy equation to solve. 

All you have to do is to recall that the derivative of an exponential is an exponential. de f /dt = e . 
Substitute this exponential for x(t), and of course it can't work as a solution; it doesn't even make 
sense dimensionally. What is e to the power of a day? You need something in the exponent to make it 
dimensionless, e at . Also, the function x is supposed to give you a position, with dimensions of length. 
Use another constant: x(t) = Ae at . Plug this into the differential equation (4.4) to find 

mAa 2 e at + bAae at + kAe at = Ae at [ma 2 + ba + k]=0 

The product of factors is zero, and the only way that a product of two numbers can be zero is if one of 
the numbers is zero. The exponential never vanishes, and for a non-trivial solution A ^ 0, so all that's 
left is the polynomial in a. 

9 , , , — b ± Vb 2 — 4km , . 

ma + oa + k = 0, with solutions a = (4.5) 

2m 



The position function is then 



x 



(t) = Ae ait + Be a2t (4.6) 



where A and B are arbitrary constants and a\ and «2 are the two roots. 

Isn't this supposed to be oscillating? It is a harmonic oscillator after all, but the exponentials 
don't look very oscillatory. If you have a mass on the end of a spring and the entire system is immersed 
in honey, it won't do much oscillating! Translated into mathematics, this says that if the constant b is 
too large, there is no oscillation. In the equation for a, if b is large enough the argument of the square 
root is positive, and both a's are real — no oscillation. Only if b is small enough does the argument of 
the square root become negative; then you get complex values for the a's and hence oscillations. 
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Push this to the extreme case where the damping vanishes: b = 0. Then a.\ = iy/k/m and 
a.i = —iy/k/m. Denote ujq = y/k/m. 

x(t) = Ae iu)ot + Be~ iuJot (4.7) 

You can write this in other forms using sines and cosines, see problem 4.10. To determine the arbitrary 
constant A and B you need two equations. They come from some additional information about the 
problem, typically some initial conditions. Take a specific example in which you start from the origin 
with a kick, x(0) = and x(0) = Vq. 

x(0) = = A + B, x(0) = v = iui A - ico B 

Solve for A and B to get A = —B = Vo/(2iuo). Then 

x(t) = -^- \e iUot - e~ iuJot ] = ^ sin w t 

As a check on the algebra, use the first term in the power series expansion of the sine function to see 
how x behaves for small t. The sine factor is sinwo^ ~ Wot, and then x(t) is approximately Vot, just as 
it should be. Also notice that despite all the complex numbers, the final answer is real. This is another 
check on the algebra. 

Damped Oscillator 

If there is damping, but not too much, then the a's have an imaginary part and a negative real part. 
(Is it important whether it's negative or not?) 

-b±i\fAkm - b 2 b . , , . / k b 2 

a = = ± tu! , where u) = \ x (4.8) 

2m 2m V m 4m z 

This represents a damped oscillation and has frequency a bit lower than the one in the undamped case. 
Use the same initial conditions as above and you will get similar results (let 7 = b/2m) 

x(t) = Ae^ +iu) '^ + Be^- iul 'V 
x(0) = A + B =0, v x (0) = (-j + iu)')A+(-j-iLj')B = v (4.9) 

The two equations for the unknowns A and B imply B = —A and 

2iu'A = v , so x(t) = —e'^ \e iuj,t - e^} = ^e"^ sinw't (4.10) 



For small values of t, the first terms in the power series expansion of this result are 

x(t) = — [1 - jt + 7 2 t 2 / 2 - • • •] [u't - u'H 3 /6 + ...}= v t - v -ft 2 + ... 

CO 

The first term is what you should expect, as the initial velocity is v x = v$. The negative sign in the 
next term says that it doesn't move as far as it would without the damping, but analyze it further. 
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Does it have the right size as well as the right sign? It is —v^t 2 = — Vo(b/2m)t 2 . But that's an 
acceleration: a x t 2 /2. It says that the acceleration just after the motion starts is a x = —bvo/m. Is 
that what you should expect? As the motion starts, the mass hasn't gone very far so the spring doesn't 
yet exert much force. The viscous friction is however —bv x . Set that equal to ma x and you see that 
— Vo'ft 2 has precisely the right value: 

9 bo 1 — bvo 9 

x(t) « v t - vojr = v t - v —t z = v t + r 

2m 2 m 

The last term says that the acceleration starts as a x = —bvo/m, as required. 

In Eq. (4.8) I assumed that the two roots of the quadratic, the two a's, are different. What if 
they aren't? Then you have just one value of a to use in defining the solution e at in Eq. (4.9). You 
now have just one arbitrary constant with which to match two initial conditions. You're stuck. See 
problem 4.11 to understand how to handle this case (critical damping). It's really a special case of what 
I've already done. 

What is the energy for this damped oscillator? The kinetic energy is mv 2 /2 and the potential 
energy for the spring is kx 2 /2. Is the sum constant? No. 

If F x = ma x = -kx + F Xifrict , then 

"Force times velocity" is a common expression for power, and this says that the total energy is decreasing 
according to this formula. For the wet friction used here, this is dE/dt = —bv 2 , and the energy 
decreases exponentially on average. 

4.2 Forced Oscillations 

What happens if the equation is inhomogeneous? That is, what if there is a term that doesn't involve 
x or its derivatives at all. In this harmonic oscillator example, apply an extra external force. Maybe it's 
a constant; maybe it's an oscillating force; it can be anything you want not involving x. 

m^ = -kx-b^ + F ext (t) (4.12) 

The key result that you need for this class of equations is very simple to state and not too difficult to 
implement. It is a procedure for attacking any linear inhomogeneous differential equation and consists 
of three steps. 

1. Temporarily throw out the inhomogeneous term [here F ext (t)] and completely solve the 
resulting homogeneous equation. In the current case that's what you just saw when I 
worked out the solution to the differential equation md 2 x/dt 2 + bdx/dt + kx = 0. 

2. Find any one solution to the full inhomogeneous equation. Note that for step one you 
have to have all the arbitrary constants present; for step two you do not. [X] n h(t)] 

3. Add the results of steps one and two. [Xhom(^) +^inh(^)] 

I've already done step one. To carry out the next step I'll start with a particular case of the 
forcing function. If F ext (t) is simple enough, you should be able to guess the answer to step two. If it's 
a constant, then a constant will work for x. If it's a sine or cosine, then you can guess that a sine or 
cosine or a combination of the two should work. If it's an exponential, then guess an exponential — 
remember that the derivative of an exponential is an exponential. If it's the sum of two terms, such 



4 — Differential Equations 



as a constant and an exponential, it's easy to verify that you add the results that you get for the two 
cases separately. If the forcing function is too complicated for you to guess a solution then there's a 
general method using Green's functions that I'll get to in section 4.6. 
Choose a specific example 

F«t(*) = F [l-e-#] (4.13) 

This starts at zero and builds up to a final value of Fq. It does it slowly or quickly depending on /3. 

F 




Start with the first term, Fq, for external force in Eq. (4.12). Try x(t) = C and plug into that 
equation to find 

kC = F 

This is simple and determines C. 

Next, use the second term as the forcing function, —Foe~P t . Guess a solution x(t) = C 'e~^ f 
and plug in. The exponential cancels, leaving 



mC/3 2 - bC'P + kC 



or 



a 



-Fa 



m/3 2 -b/3 + k 

The total solution for the inhomogeneous part of the equation is then the sum of these two expressions. 

1 1 



X\nh(t) = F 



-pt 



k m(3 2 -b/3 + k 

The homogeneous part of Eq. (4.12) has the solution found in Eq. (4.6) and the total is 



x(t) = x hom (t) + x- mh (t) = X 



(t) = Ae ait + Be° 2t + F Q 



1 



-6t 



k m(3 2 -b(3 + k 



(4.14) 



There are two arbitrary constants here, and this is what you need because you have to be able to 
specify the initial position and the initial velocity independently; this is a second order differential 
equation after all. Take for example the conditions that the initial position is zero and the initial 
velocity is zero. Everything is at rest until you start applying the external force. This provides two 
equations for the two unknowns. 



x(0) = = A + B + F 



m(3 2 — b(3 



x(0) = = Aai + Ba 2 + F 



k(m(3 2 -b(3 + k) 

(3 



m(3 2 -b/3 + k 



Now all you have to do is solve the two equations in the two unknowns A and B. Take the first, 
multiply it by a% and subtract the second. This gives A. Do the same with at\ instead of 0(2 to get B. 
The results are 



.4 



1 ^a 2 (m{3 2 -b(3)-k(3 



OL\ — «2 



k(m(3 2 -bp + k) 
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Interchange ot\ and a 2 to get B. 
The final result is 



F (a 2 (m/3 2 - b(5) - fc/5) e Ql * - (ax {m(3 2 - b(5) - fc/3) e« 2 * 



+ Fo[ l T - ^ m \ a , (1.15) 



«i — «2 k(m(3 2 — b[3 + k) 

'1 1 

A; m(3 2 -b(3 + k 

If you think this is messy and complicated, you haven't seen messy and complicated. When it takes 
20 pages to write out the equation, then you're entitled say that it is starting to become involved. 

Why not start with a simpler example, one without all the terms? The reason is that a complex 
expression is often easier to analyze than a simple one. There are more things that you can do to it, and 
so more opportunities for it to go wrong. The problem isn't finished until you've analyzed the supposed 
solution. After all, I may have made some errors in algebra along the way. Also, analyzing the solution 
is the way you learn how these functions work. 

1. Everything in the solution is proportional to Fq and that's not surprising. 

2. I'll leave it as an exercise to check the dimensions. 

3. A key parameter to vary is (3. What should happen if it is either very large or very 
small? In the former case the exponential function in the force drops to zero quickly 
so the force jumps from zero to Fq in a very short time — a step in the limit that 

/3->0. 

4. If /3 is very small the force turns on very gradually and gently, as though you are being 
very careful not to disturb the system. 

Take point 3 above: for large /3 the dominant terms in both numerator and denominator every- 
where are the m/3 2 terms. This result is then very nearly 

, t) ^ Fq (« 2 (m/3 2 ))e^-(a 1 (m/3 2 ))e^ + / 1 _ _J_ g _ /W 
«i — ct2 km(3 2 \k (m(3 2 ) 



Ft 



o 



[(a 2 e ait - axe"*] + F Q \ 



k(ax — ot-2) 
Use the notation of Eq. (4.9) and you have 

x(t) w — —, ^-v [((-7 - iu')e^ +iuj ^ - (-7 + iu')e^- lbJ 'n + F \- 

w fc(- 7 + «£*;'- (-7-zu/)) LV ; J k 

= ° [ - ^ sin u/t - cos co't] + F T (4.16) 

K CO rC 
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At time t = this is still zero even with the approximations. That's comforting, but if it hadn't 
happened it's not an insurmountable disaster. This is an approximation to the exact answer after all, 
so it could happen that the initial conditions are obeyed only approximately. The exponential terms 
have oscillations and damping, so the mass oscillates about its eventual equilibrium position and after 
a long enough time the oscillations die out and you are left with the equilibrium solution x = F^/k. 

Look at point 4 above: For small f3 the /3 2 terms in Eq. (4.15) are small compared to the (3 terms 
to which they are added or subtracted. The numerators of the terms with e at are then proportional to 
(3. The denominator of the same terms has a k — b/3 in it. That means that as f3 — > 0, the numerator 
of the homogeneous term approaches zero and its denominator doesn't. The last terms, that came 
from the inhomogeneous part, don't have any f3 in the numerator so they don't vanish in this limit. 
The approximate final result then comes solely from the X\ n h(t) term. 

It doesn't oscillate at all and just gradually moves from equilibrium to equilibrium as time goes on. It's 
what you get if you go back to the differential equation (4.12) and say that the acceleration and the 
velocity are negligible. 

dp x dx 1 

m-Tp[nO] = -kx - b^[^ 0] + F ext (t) =^ x^j_F ext (t) 

The spring force nearly balances the external force at all times; this is "quasi-static," in which the 
external force is turned on so slowly that it doesn't cause any oscillations. 

4.3 Series Solutions 

A linear, second order differential equation can always be rearranged into the form 

y" + P{x)y' + Q(x)y = R(x) (4.17) 

If at some point Xq the functions P and Q are well-behaved, if they have convergent power series 
expansions about Xq, then this point is called a "regular point" and you can expect good behavior of 
the solutions there — at least if R is also regular there. 

I'll look just at the case for which the inhomogeneous term R = 0. If P or Q has a singularity 
at Xq, perhaps something such as l/(x — Xo) or yjx — Xq, then Xq is called a "singular point" of the 
differential equation. 

Regular Singular Points 

The most important special case of a singular point is the " regular singular point" for which the behaviors 
of P and Q are not too bad. Specifically this requires that {x — Xq)P(x) and (x — Xq) 2 Q(x) have no 
singularity at Xq. For example 

y " + l y, + ^ y = and y " + ^y' +x y = 

have singular points at X = 0, but the first one is a regular singular point and the second one is not. 
The importance of a regular singular point is that there is a procedure guaranteed to find a solution near 
a regular singular point (Frobenius series). For the more general singular point there is no guaranteed 
procedure (though there are a few tricks* that sometimes work). 

* The book by Bender and Orszag: "Advanced mathematical methods for scientists and engineers" 
is a very readable source for this and many other topics. 
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Examples of equations that show up in physics problems are 

y" + y = o 

(1 - x 2 )y" - 2xy' + t{l + l)y = 
x 2 y" + xy' + (x 2 - n 2 )y = 
xy" + (a + 1 — x)y' + ny = 

These are respectively the classical simple harmonic oscillator, Legendre equation, Bessel equation, 
generalized Laguerre equation. 

A standard procedure to solve these equations is to use series solutions, but not just the standard 
power series such as those in Eq. (2.4). Essentially, you assume that there is a solution in the form of 
an infinite series and you systematically compute the terms of the series. I'll pick the Bessel equation 
from the above examples, as the other three equations are done the same way. The parameter n in that 
equation is often an integer, but it can be anything. It's common for it to be 1/2 or 3 /2 or sometimes 
even imaginary, but there's no need to make any assumptions about it for now. 

Assume a solution in the form : 



regular singular points at ±1 

(4. 18) 

regular singular point at zero 
regular singular point at zero 



00 



Frobenius Series: y(x) = ^ a k x k+s (a ^0) (4.19) 



If s = or a positive integer, this is just the standard Taylor series you saw so much of in chapter two, 
but this simple-looking extension makes it much more flexible and suited for differential equations. It 
often happens that s is a fraction or negative, but this case is no harder to handle than the Taylor 
series. For example, what is the series expansion of (cosx)/:r about the origin? This is singular at 
zero, but it's easy to write the answer anyway because you already know the series for the cosine. 

cos a; 1 x x 3 x 5 

+ ;tt - zr- + ••• 



x a; 2 24 720 

It starts with the term l/x corresponding to s = — 1 in the Frobenius series. 

Always assume that ^ 0, because that just defines the coefficient of the most negative power, 
X s . If you allow it be zero, that's just the same as redefining s and it gains nothing except confusion. 
Plug this into the Bessel differential equation. 

x 2 y" + xy' + {x 2 -n 2 )y = 



00 00 00 

x 2 ^2 a k(k + s)(k + s - l)x k+s ~ 2 + x ^2 a k(k + s)x k+s ' 1 + (x 2 - n 2 ) ^ a k x k+s = 

fc=0 fc=0 k=Q 

00 00 00 00 

a k(k + s)(k + s- l)x k+s + ^2a k (k + s)x k+s + ^ a k x k+s+2 - n 2 ^ a k x k+s = 

k=0 k=0 k=0 

00 00 

J2 a k [(k + s)(k + s - 1) + (k + s) - n 2 ]x k+s + ^ a k x k+s+2 = 

k=0 k=0 



The coefficients of all the like powers of x must match, and in order to work out the matches efficiently, 
and so as not to get myself confused in a mess of indices, I'll make an explicit change of the index in 
the sums. Do this trick every time. It keeps you out of trouble. 
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Let £ = k in the first sum. Let £ = k + 2 in the second. Explicitly show the limits of the index 
on the sums, or you're bound to get it wrong. 

oo oo 

H [{£ + s) 2 - n 2 ] x l+s + £ a e _ 2 x e+s = 

£=0 i=2 

The lowest power of x in this equation comes from the £ = term in the first sum. That coefficient of 
X s must vanish, (do ^ 0) 

a [s 2 -n 2 ] = (4.20) 

This is called the indicia! equation. It determines s, or in this case, maybe two s's. After this, set to 
zero the coefficient of x^ +s . 

a £ [(£ + s) 2 - n 2 ] + a e _ 2 = (4.21) 
This determines a 2 in terms of <2o; it determines in terms of a 2 etc. 



"' ~"'--'0 + s )2 _ n V ^~ 2 ' 4 ' ••• 



For example, if n = 0, the indicial equation says s = 0. 

1 1 1 1 



22' "^42 |U/U 2 2 4 2 ' 6 2 u 2 2 4 2 6 2 



1 °° (x/2) 2k 

a 2 k = (-l) fc a ^2fcfcj2 then l/(^) = «o^(" 1 ) ( fch2 = Q o^o(x) (4.22) 

fc=o ^ 



and in the last equation I rearranged the factors and used the standard notation for the Bessel function, 

J n {x). 

This is a second order differential equation. What about the other solution? This Frobenius 
series method is guaranteed to find one solution near a regular singular point. Sometimes it gives both 
but not always, and in this example it produces only one. There are procedures that will let you find 
the second solution to this sort of second order differential equation. See problem 4.49 for one such 
method. 

For the case n = 1 fe the calculations just above will produce two solutions. The indicial equation 
gives s = i 1 /^- After that, the recursion relation for the coefficients give 

1 1 1 1 

a £ — - a l-2,f, , „no —o — ~ a £-2 ffo , n /,„ — - a £-2 fU (, , „\ _ _C ^-2" 



sf-n 2 ^-H 2 + 2£s l ~ z £(£ + 2s) l ~ z £(£±l) 
For the s = +V2 result 

1 1 1 

0-2 = ~Q<o 7. — 7; «4 = -a-2-r- r = +«o 



(2k + iy. 

This solution is then 



2-3 ' 4-5 "2-3-4-5 

&2k = (-l) fc ao- 



sk, 1 



y(x) = aox 1 / 2 



2 A 
X X 



1 -¥ + 5! 
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This series looks suspiciously like the series for the sine function, but is has some of the x's or some of 
the factorials in the wrong place. You can fix that if you multiply the series in brackets by x. You then 
have 

y(x) = aox^ 1 / 2 



X ~3l + 5! 



sm x , , s 



I'll leave it to problem 4.15 for you to find the other solution. 

Do you need to use a Frobenius series instead of just a power series for all differential equations? 
No, but I recommend it. If you are expanding about a regular point of the equation then a power series 
will work, but I find it more systematic to use the same method for all cases. It's less prone to error. 

4.4 Some General Methods 

It is important to be familiar with the arsenal of special methods that work on special types of differential 
equations. What if you encounter an equation that doesn't fit these special methods? There are some 
techniques that you should be familiar with, even if they are mostly not ones that you will want to 
use often. Here are a couple of methods that can get you started, and there's a much broader set of 
approaches under the heading of numerical analysis; you can explore those in section 11.5. 

If you have a first order differential equation, dx/dt = /(£,£), with initial condition x(to) = Xq 
then you can follow the spirit of the series method, computing successive orders in the expansion. 
Assume for now that the function / is smooth, with as many derivatives as you want, then use the 
chain rule a lot to get the higher derivatives of x 

d 2 x _df df dx _ .. _ f . . 
W_~dt + dxdt~ X ~ ft + fxX 

X = f tt + 2fxt± + jxxX 1 + f x 'x = fu + 2f xt x + f xx x 2 + f x [f t + f x A} 
x{t) = x Q + f(x , t )(t - t ) + lx(t )(t - t ) 2 + lx(t )x(t )(t - t ) 3 + ■■■ (4.24) 

Here the dot-notation (i etc.) is a standard shorthand for derivative with respect to time. This is 
unlike using a prime for derivative, which is with respect to anything you want. These equations show 
that once you have the initial data (to,Xo), you can compute the next derivatives from them and from 
the properties of /. Of course if / is complicated this will quickly become a mess, but even then it can 
be useful to compute the first few terms in the power series expansion of x. 
For example, x = f(x, t) = Ax 2 (l + ut) with to = and Xq = a. 

xo = Aa 2 , x = Aa 2 uj + 2A 2 a 3 , x = AA 2 a 3 u + 2A 3 a 4 + 2Aa [Aa 2 u + 2A 2 a 3 ] (4.25) 

Ifv4 = 1/ m-sand uj = l/s with a = 1 m this is 

x(t) = 1 + t + \t 2 + 2t 3 + • • • 

You can also solve this example exactly and compare the results to check the method. 

What if you have a second order differential equation? Pretty much the same thing, though it is 
sometimes convenient to make a slight change in the appearance of the equations when you do this. 

x = f(x,i,t) can be written x = v, v = f(x, v,t) (4-26) 

so that it looks like two simultaneous first order equations. Either form will let you compute the higher 
derivatives, but the second one often makes for a less cumbersome notation. You start by knowing to. 

Xo, and now vq = Xq. 
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Some of the numerical methods you will find in chapter 11 start from the ideas of these expansions, 
but then develop them along different lines. 

There is an iterative methods that of more theoretical than practical importance, but it's easy 
to understand. I'll write it for a first order equation, but you can rewrite it for the second (or higher) 
order case by doing the same thing as in Eq. (4.26). 

x = f(x,t) with a; (to) = generates Xi(t) = / dt' f(xo,t') 

J to 

This is not a solution of the differential equation, but it forms the starting point to find one because 
you can iterate this approximate solution X\ to form an improved approximation. 

x k (t)= f t dff(x k _ l (1/),f), k = 2,3,... (4.27) 

This will form a sequence that is usually different from that of the power series approach, though the 
end result better be the same. This iterative approach is used in one proof that shows under just what 
circumstances this differential equation x = f has a unique solution. 

4.5 Trigonometry via ODE's 

The differential equation u" = —u has two independent solutions. The point of this exercise is to derive 
all (or at least some) of the standard relationships for sines and cosines strictly from the differential 
equation. The reasons for spending some time on this are twofold. First, it's neat. Second, you have 
to get used to manipulating a differential equation in order to find properties of its solutions. This is 
essential in the study of Fourier series as you will see in section 5.3. 

Two solutions can be defined when you specify boundary conditions. Call the functions c(x) and 
s(x), and specify their respective boundary conditions to be 

c(0) = 1, c'(0) = 0, and s(0) = 0, s'(0) = 1 (4.28) 

What is s'(x)7 First observe that s' satisfies the same differential equation as s and c: 

u" = —u => (u 1 )" = (u")' = —vf, and that shows the desired result. 

This in turn implies that s' is a linear combination of s and c, as that is the most general solution to 
the original differential equation. 

s'(x) = Ac(x) +Bs(x) 

Use the boundary conditions: 

s '(0) = 1 = Ac(0) + Bs{0) = A 
From the differential equation you also have 

s"(0) = -s(0) = = Ac'(0) + Bs' (0) = B 

Put these together and you have 

s'(x) = c(x) And a similar calculation shows c'(x) = —s(x) (4.29) 

What is c(x) 2 + s(x) 2 l Differentiate this expression to get 

^[c(a;) 2 + s(x) 2 } = 2c(x)c'(x) + 2s(x)s'(x) = -2c(x)s(x) + 2s(x)c(x) = 

This combination is therefore a constant. What constant? Just evaluate it at x = and you see that 
the constant is one. There are many more such results that you can derive, but that's left for the 
exercises, problem 4.21 et seq. 
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4.6 Green's Functions 

Is there a general way to find the solution to the whole harmonic oscillator inhomogeneous differential 
equation? One that does not require guessing the form of the solution and applying initial conditions? 
Yes there is. It's called the method of Green's functions. The idea behind it is that you can think of 
any force as a sequence of short, small kicks. In fact, because of the atomic nature of matter, that's 
not so far from the truth. If you can figure out the result of an impact by one molecule, you can add 
the results of many such kicks to get the answer for 10 23 molecules. 

I'll start with the simpler case where there's no damping, b = in the harmonic oscillator 
equation. 

mx + kx = F ext (t) (4.30) 

Suppose that everything is at rest at the origin and then at time if the external force provides a small 
impulse. The motion from that point on will be a sine function starting at t', 

A sin (w (t-tf)) (t>f) (4.31) 

The amplitude will depend on the strength of the kick. A constant force F applied for a very short 
time, At', will change the momentum of the mass by mAv x = FAt'. If this time interval is short 
enough the mass doesn't have a chance to move very far before the force is turned off, then from that 
time on it's subject only to the — kx force. This kick gives m a velocity FAt' /m, and that's what 
determines the unknown constant A. 

Just after t = t' , v x = Auq = FAt'/m. This determines A, so the position of m is 

FA % sin {ou (t-t')) (t>t>) 
(t < t>) 



impulse 



x(t) 



mujo 




(4.32) 




When the external force is the sum of two terms, the total solution is the sum of the solutions 
for the individual forces. If an impulse at one time gives a solution Eq. (4.32), an impulse at a later 
time gives a solution that starts its motion at that later time. The key fact about the equation that 
you're trying to solve is that it is linear, so you can get the solution for two impulses simply by adding 
the two simpler solutions. 

d 2 Xi 



ni- 



di 2 



+ kx\ = Fi(t) and 



d 2 x 2 , , 



m) 



then 



m 



d 2 (xi +x 2 ) 



k{x 1 + x 2 ) = F 1 (t) + F 2 (t) 



vwww 



X\ + x 2 



+ 



X 2 \f\ 






V 




/v\ 


l\ 



Vvwwv 
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The way to make use of this picture is to take a sequence of contiguous steps. One step follows 
immediately after the preceding one. If two such impulses are two steps 



F{t ) (to<*<ti) 
(elsewhere) 



and F 1 = l™ (*}<«<« 
[ (elsewhere) 



mi + kx = F + Fx 



(4.33) 



then if Xo is the solution to Eq. (4.30) with only the Fq on its right, and X\ is the solution with only 
F\, then the full solution to Eq. (4.33) is the sum, Xq + X\. 

Think of a general forcing function F x>ext (t) in the way that you would set up an integral. 
Approximate it as a sequence of very short steps as in the picture. Between t k and t k+1 the force is 
essentially F{t k ). The response of m to this piece of the total force is then Eq. (4.32). 



x k (t) 



^MMk s in(co (t-t k )) (t>t k ) 
3 (t < t k ) 




+ 



+ 



To complete this idea, the external force is the sum of a lot of terms, the force between t\ and 
ti, that between ti and t-$ etc. The total response is the sum of all these individual responses. 



x(t) = yi^ ks[n (^(t-t k )) (t>t k ) 



k 



For a specified time t, only the times t k before and up to t contribute to this sum. The impulses 
occurring at the times after the time t can't change the value of x(t); they haven't happened yet. In 
the limit that At k — > 0, this sum becomes an integral. 



x(t) 



dt ,F ^ sin (wo(*-<0) 
-oo m ^o 



(4.34) 



Apply this to an example. The simplest is to start at rest and begin applying a constant force 
from time zero on. 



F ext (t) 



F (t>0) 
(t < 0) 



x(t) 



J-t-l . 
at sin 



M*-*0) 



and the last expression applies only for t > 0. It is 

F 



x{t) 



[l — cos(cuot)] 



(4.35) 
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As a check for the plausibility of this result, look at the special case of small times. Use the power 
series expansion of the cosine, keeping a couple of terms, to get 

x(t) - A [1 " (1 " ("ot?/2)] = Ai" 2 = -f 
mwg L v /J mu^ 2 m 2 

and this is just the result you'd get for constant acceleration Fo/m. In this short time, the position 
hasn't changed much from zero, so the spring hasn't had a chance to stretch very far, so it can't apply 
much force, and you have nearly constant acceleration. 

This is a sufficiently important subject that it will be repeated elsewhere in this text. A completely 
different approach to Green's functions will appear is in section 15.5, and chapter 17 is largely devoted 
to the subject. 

4.7 Separation of Variables 

If you have a first order differential equation — I'll be more specific for an example, in terms of x and 
t — and if you are able to move the variables around until everything involving x and dx is on one 
side of the equation and everything involving t and dt is on the other side, then you have "separated 
variables." Now all you have to do is integrate. 

For example, the total energy in the undamped harmonic oscillator is E = mv 2 /2 + kx 2 /2. 
Solve for dx/dt and 

Tt = ^ E - fa2 /2) <" 6 > 

To separate variables, multiply by dt and divide by the right-hand side. 

= = dt 



j2_ 
m 1 



(E-kx 2 /2) 

Now it's just manipulation to put this into a convenient form to integrate. 

fm dx ^ f dx f fk ^ 

Make the substitution x = asm.9 and you see that if a 2 = 2E/k then the integral on the left simplifies. 

acos9d9 f k ,, „ . x 



f acosOdO f k „ „ _■, x 

/ — . = / \ — dt so v = sin — = coot + C 

J aJl-sm 2 9 J V m a 



Vl - sin 2 6 

or x(t) = asm(cuot + C) where tu = \jkjm 

An electric circuit with an inductor, a resistor, and a battery has a differential equation for the 
current flow: 

L^ + IR = V (4.37) 

1 I 

Manipulate this into 

L% = V -IR, then L v ^ m =dt 
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Now integrate this to get 

L J VQ dI IR =t + C, or - ^]n(Vo- IR)=t + C 
Solve for the current / to get 

RI(t) =V - e-iL/m+C) ( 438 ) 

Now does this make sense? Look at the dimensions and you see that it doesn't, at least not yet. The 
problem is the logarithm on the preceding line where you see that its units don't make sense either. 
How can this be? The differential equation that you started with is correct, so how did the units get 
messed up? It goes back to the standard equation for integration, 



J dx/x = In x + C 



If x is a length for example, then the left side is dimensionless, but this right side is the logarithm of 
a length. It's a peculiarity of the logarithm that leads to this anomaly. You can write the constant of 
integration as C = — lnC" where C is another arbitrary constant, then 



/x 
dx/x = lax + C = Inx — lnC = In 



If C is a length this is perfectly sensible dimensionally. To see that the dimensions in Eq. (4.38) will 
work themselves out (this time), put on some initial conditions. Set 1(0) = so that the circuit starts 
with zero current. 

R.0 = V - e-(W0+C) impNes e -(L/R)(C) = y Q 

RI(t) = V - V e~ Lt / R or I(t) = (1 - e~ Lt / R )V /R 

and somehow the units have worked themselves out. Logarithms do this, but you still better check. 
The current in the circuit starts at zero and climbs gradually to its final value I = Vq/R. 

4.8 Circuits 

The methods of section 4.1 apply to simple linear circuits, and the use of complex algebra as in that 
section leads to powerful and simple ways to manipulate such circuit equations. You probably remember 
the result of putting two resistors in series or in parallel, but what about combinations of capacitors or 
inductors under the same circumstances? And what if you have some of each? With the right tools, 
all of these questions become the same question, so it's not several different techniques, but one. 

If you have an oscillating voltage source (a wall plug), and you apply it to a resistor or to a 
capacitor or to an inductor, what happens? In the first case, V = IR of course, but what about the 
others? The voltage equation for a capacitor is V = q/C and for an inductor it is V = Ldl/dt. A 
voltage that oscillates at frequency u is V = Vocosuit, but using this trigonometric function forgoes 
all the advantages that complex exponentials provide. Instead, assume that your voltage source is 
V = Voe lbjt with the real part understood. Carry this exponential through the calculation, and take 
the real part only at the end — often you won't even need to do that. 



— o 

— T/ ijr^SffiSS& 



-0 
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These are respectively 



V e iu3t = IR = I e iLut R 

V e iujt = q/C =► iuV e iut = q/C = I/C = I e lujt /C 
V e iujt = LI = iuLI = iuLI e iu}t 

In each case the exponential factor is in common, and you can cancel it. These equations are then 

V = IR V = I/zuC V = icoLI 

All three of these have the same form: V = (something times)/, and in each case the size of the 
current is proportional to the applied voltage. The factors of i implies that in the second and third 
cases the current is ±90° out of phase with the voltage cycle. 

The coefficients in these equations generalize the concept of resistance, and they are called 
"impedance," respectively resistive impedance, capacitive impedance, and inductive impedance. 



V = Z R I = RI 



V = Z C I 



1 



ituC 



V = Z L I = ujLI 



(4.39) 



Impedance appears in the same place as does resistance in the direct current situation, and this implies 
that it can be manipulated in the same way. The left figure shows two impedances in series. 







z 2 









r Z 1 -4 



L- Z, -I 



The total voltage from left to right in the left picture is 

V = ZJ + Z 2 I = (Z 1 + Z 2 )I = Z totai I 



(4.40) 



It doesn't matter if what's inside the box is a resistor or some more complicated impedance, it matters 
only that each box obeys V = ZI and that the total voltage from left to right is the sum of the two 
voltages. Impedances in series add. You don't need the common factor e %Lut . 

For the second picture, for which the components are in parallel, the voltage is the same on each 
impedance and charge is conserved, so the current entering the circuit obeys 



/ = h + I 2 , 



then 



V 



J total 



— + — 

Z\ z 2 



or 



J total 



1 1 

Z\ Z 2 



(4.41) 



Impedances in parallel add as reciprocals, so both of these formulas generalize the common equations 
for resistors in series and parallel. They also include as a special case the formula you may have seen 
before for adding capacitors in series and parallel. 

In the example Eq. (4.37), if you replace the constant voltage by an oscillating voltage, you have 
two impedances in series. 



Z R + Z L = R + icoL => I = V / (R + iuL) 
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What happened to the e~ Lt l R term of the previous solution? This impedance manipulation tells you 
the inhomogeneous solution; you still must solve the homogeneous part of the differential equation and 
add that. 

L ^L + ir = o l{t) = Ae- Rt l L 



dt 



The total solution is the sum 



I(t) 
real part 



Ae -Rt/L + yr 



OP 



R + io;L 
VR 2 +uj 2 L 2 



where 



tan 



_i 

IT 



(4.42) 



How did that last manipulation come about? Change the complex number R + iuL in the denominator 
from rectangular to polar form. Then the division of the complex numbers becomes easy. The dying 
exponential is called the "transient" term, and the other term is the "steady-state" term. 
The denominator is 



R + iuL = a + i(3 = a 2 + /3 2 



a + i(3 



(4.43) 



The reason for this multiplication and division by the same factor is 
that it makes the final fraction have magnitude one. That allows me 
to write it as an exponential, e 1 ^ . From the picture, the cosine and the 
sine of the angle <p are the two terms in the fraction. 











a 



a + i/3 = \J a 2 + /3 2 (cos + 



f3 2 e** 



and 



tan (f) = /3/a 



In summary, 



V = IZ 



V 

z 



\Z\e i(t> 



V 



VR 2 + u 2 L 2 e i( t ) 



To satisfy initial conditions, you need the parameter A, but you also see that it gives a dying 
exponential. After some time this transient term will be negligible, and only the oscillating steady-state 
term is left. That is what this impedance idea provides. 

In even the simplest circuits such as these, that fact that Z is complex implies that the applied 
voltage is out of phase with the current. Z = \Z\e 1 ^, so I = V/Z has a phase change of —(f) from V. 

What if you have more than one voltage source, perhaps the second having a different frequency 
from the first? Remember that you're just solving an inhomogeneous differential equation, and you are 
using the methods of section 4.2. If the external force in Eq. (4.12) has two terms, you can handle 
them separately then add the results. 



4.9 Simultaneous Equations 

What's this doing in a chapter on differential equations? Patience, 
unknowns: 



Solve two equations in two 



dx(X) - 6x(Y): 



(X) ax + by = e 
(Y) cx + dy = f 

Similarly, multiply (Y) by a and (X) by c and subtract: 

acx + ady — acx — cby 
(ad — bc)y 



adx + bdy 



- box — bdy 
(ad — bc)x 



ed 
ed 



fb 
fb 



fa 
fa 



ec 
ec 
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x = 



ed — fb 
ad — be 



V = 



fa — ec 



ad — be 



(4.44) 



provided that ad — be ^ 0. This expression appearing in both denominators is the determinant of the 
equations. 

Classify all the essentially different cases that can occur with this simple-looking set of equations 
and draw graphs to illustrate them. If this looks like problem 1.23, it should. 



1. The solution is just as in Eq. (4.44) above and nothing goes wrong. There is exactly one 
solution. The two graphs of the two equations are two intersecting straight lines. 

2. The denominator, the determinant, is zero and the numerator isn't. This is impossible and 
there are no solutions. When the determinant vanishes, the two straight lines are parallel and the fact 
that the numerator isn't zero implies that the two lines are distinct and never intersect. (This could also 
happen if in one of the equations, say (X), a = b = and e ^ 0. For example = 1. This obviously 
makes no sense.) 

3a. The determinant is zero and so are both numerators. In this case the two lines are not only 
parallel, they are the same line. The two equations are not really independent and you have an infinite 
number of solutions. 

3b. You can get zero over zero another way. Both equations (X) and (Y) are = 0. This sounds 
trivial, but it can really happen. Every x and y will satisfy the equation. 

4. Not strictly a different case, but sufficiently important to discuss it separately: suppose that 
the right-hand sides of (X) and (Y) are zero, e = / = 0. If the determinant is non-zero, there is a 
unique solution and it is x = 0, y = 0. 

5. With e = / = 0, if the determinant is zero, the two equations are the same equation and 
there are an infinite number of non-zero solutions. 

In the important case for which e = / = and the determinant is zero, there are two cases: 
(3b) and (5). In the latter case there is a one-parameter family of solutions and in the former case 
there is a two-parameter family. Put another way, for case (5) the set of all solutions is a straight line, 
a one-dimensional set. For case (3b) the set of all solutions is the whole plane, a two-dimensional set. 





y 



X 



3b. 



Example: consider the two equations 



kx + (k- l)y = 0, (1 - k)x + {k - ify = 
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For whatever reason, I would like to get a non-zero solution for x and y. Can I? The condition depends 
on the determinant, so take the determinant and set it equal to zero. 

k(k - I) 2 - {1 - k)(k - 1) = 0, or (A; + 1)(A;-1) 2 = 

There are two roots, k = — 1 and k = +1. In the k = — 1 case the two equations become 

-x - 2y = 0, and 2x + 4y = 

The second is just —2 times the first, so it isn't a separate equation. The family of solutions is all those 
x and y satisfying x = —2y, a straight line. 
In the k = +1 case you have 

x + 0y = 0, and = 
The solution to this is x = and y = anything and it is again a straight line (the y-axis). 
4.10 Simultaneous ODE's 

Single point masses are an idealization that has some application to the real world, but there are many 
more cases for which you need to consider the interactions among many masses. To approach this, 
take the first step, from one mass to two masses. 

Two masses are connected to a set of springs and fastened be- ^ k k 

tween two rigid walls as shown. The coordinates for the two masses ■ 1 3 2 ■ 

(moving along a straight line) are X\ and x 2 , and I'll pick the zero p^&H — 
point for these coordinates to be the positions at which everything is at L^. L^. 

equilibrium — no total force on either. When a mass moves away from Xl x<2 

its equilibrium position there is a force on it. On m\, the two forces are 

proportional to the distance by which the two springs k\ and k 3 are stretched. These two distances are 
X\ and X\ — x 2 respectively, so F x = ma x applied to each mass gives the equations 

mi^^- = -hxi - k 3 (x! - x 2 ), and m 2 ^ = -k 2 x 2 - k 3 (x 2 - 21) (4.45) 

I'm neglecting friction simply to keep the algebra down. These are linear, constant coefficient, homo- 
geneous equations, just the same sort as Eq. (4.4) except that there are two of them. What made 
the solution of (4.4) easy is that the derivative of an exponential is an exponential, so that when you 
substituted x(t) = Ae at all that you were left with was an algebraic factor — a quadratic equation in 
a. Exactly the same method works here. 

The only way to find out if this is true is to try it. The big difference is that there are two 
unknowns instead of one, and the amplitude of the two motions will probably not be the same. If one 
mass is a lot bigger than the other, you expect it to move less. 

Try the solution 

Xl (t) = Ae a \ x 2 {t) = Be at (4.46) 

When you plug this into the differential equations for the masses, all the factors of e at cancel, just the 
way it happens in the one variable case. 

m x o?A = -kiA - k 3 (A - B), and m 2 a 2 B = -k 2 B - k 3 (B - A) (4.47) 

Rearrange these to put them into a neater form. 

(k l + k 3 + m l a 2 )A+ (-k 3 )B = 

(-k 3 )A+(k 2 + k 3 + m 2 a 2 )B = (4.48) 
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The results of problem 1.23 and of section 4.9 tell you all about such equations. In particular, 
for the pair of equations ax + by = and cx + dy = 0, the only way to have a non-zero solution for 
x and y is for the determinant of the coefficients to be zero: ad — be = 0. Apply this result to the 
problem at hand. Either A = and B = with a trivial solution or the determinant is zero. 

(A;i + h + mi« 2 ) (k 2 + h + m 2 a 2 ) - (k 3 ) 2 = (4.49) 

This is a quadratic equation for a 2 , and it determines the frequencies of the oscillation. Note the plural 
in the word frequencies. 

Equation (4.49) is just a quadratic, but it's still messy. For a first example, try a special, 
symmetric case: rri\ = m 2 = m and k\ = k 2 . There's a lot less algebra. 

(h + k 3 + ma 2 ) 2 - {k 3 ) 2 = (4.50) 

You could use the quadratic formula on this, but why? It's already set up to be factored. 

(ki + k 3 + ma 2 - k 3 )(kx + k 3 + ma 2 + k 3 ) = 

The product is zero, so one or the other factors is zero. These determine the as. 

o ki o ki -\- 1k 3 . 

a? = — - and a\ = — — 6 4.51 

m m 

These are negative, and that's what you should expect. There's no damping and the springs provide 
restoring forces that should give oscillations. That's just what these imaginary a's provide. 

When you examine the equations ax+by = and cx+dy = the condition that the determinant 
vanishes is the condition that the two equations are really just one equation, and that the other is not 
independent of it; it is actually a multiple of the first. You must solve that equation for x and y. Here, 
arbitrarily pick the first of the equations (4.48) and find the relation between A and B. 

aj = -— (k 1 + k 3 + m(-(k 1 /m)))A+(-k 3 )B = => B = A 

a 2 = _ki + 2h ^ {k 1 + k 3 + m(-(ki + 2k 3 /m)))A+(-k 3 )B = =► B = -A 

For the first case, ot\ = = ±iy/ki/m, there are two solutions to the original differential equations. 

These are called "normal modes." 

Xl (t) = A ie iuJlt xi(t) = A 2 e~ iuit 

x 2 (t) = A ie liUlt x 2 (t) = A 2 e~ iuJlt 

The other frequency has the corresponding solutions 

xi(t) = A 3 e iU2t ^ d Xl (t) = A 4 e- iLU2t 
x 2 (t) = -A 3 e iu}2t ^ x 2 {t) = -A 4 e~ iuj2t 

The total solution to the differential equations is the sum of all four of these. 

x x (t) = A ie iuJlt + A 2 e~ i0Jlt + A 3 e lU2t + A 4 e~ luJ2t 

x 2 {t) = A ie iuJlt + A 2 e~ i0Jlt - A 3 e iuJ2t - A A e~ lW2t (4.52) 
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The two second order differential equations have four arbitrary constants in their solution. You 
can specify the initial values of two positions and of two velocities this way. As a specific example 
suppose that all initial velocities are zero and that the first mass is pushed to coordinate Xq and 
released. 

zi(0) = x = A x + A 2 + A 3 + A 4 
x 2 (o) = o = Ax + A 2 - A 3 - A A 
^xi(O) = = iuj\A\ — iuiA 2 + iu 2 A 3 — iu 2 A± 

v X 2(0) = = iu\Ai — iuj\A 2 — iu 2 A 3 + %uj 2 A^ (4.53) 
With a little thought (i.e. don't plunge blindly ahead) you can solve these easily. 

A 1 = A 2 = A 3 = A 4 = ^ 



x 2 { 



i(t) = ?j [e iuJlt + e~ ilJlt + e iLU2t + e"*" 2 *] = y [coswit + cosu 2 t] 
2(t) = °j [e iuJlt + e~ iuJlt - e iLJ2t - e~ iuJ2t ] = ^ [cos^it - cosu 2 t] 
From the results of problem 3.34, you can rewrite these as 

m /wa + WiA (uj 2 -uji 

X\{t) = XqCOS I 1 I COS I : 1 



x 2 (t) = Xosm | 1 I sin 



2 

W 2 - Ui 



(4.54) 



As usual you have to draw some graphs to understand what these imply. If the center spring k 3 
is a lot weaker than the outer ones, then Eq. (4.51) implies that the two frequencies are close to each 
other and so \ui — u 2 \ <C lO\ + cu 2 . Examine Eq. (4.54) and you see that one of the two oscillating 
factors oscillate at a much higher frequency than the other. To sketch the graph of x 2 for example you 
should draw one factor [sin {{pj 2 + Ui)t/2)] and the other factor [sin ((oj 2 — LO\)t / 2)] and graphically 
multiply them. 




The mass m 2 starts without motion and its oscillations gradually build up. Later they die down 
and build up again (though with reversed phase). Look at the other mass, governed by the equation 
for Xi(t) and you see that the low frequency oscillation from the (u 2 — lO\)/2 part is big where the one 
for x 2 is small and vice versa. The oscillation energy moves back and forth from one mass to the other. 

4.11 Legendre's Equation 

This equation and its solutions appear when you solve electric and gravitational potential problems 
in spherical coordinates [problem 9.20]. They appear when you study Gauss's method of numerical 
integration [Eq. (11.27)] and they appear when you analyze orthogonal functions [problem 6.7]. Because 
it shows up so often it is worth the time to go through the details in solving it. 



[(1 - x 2 )y'} ' + Cy = 0, or (1 - x 2 )y" - 2xy' + Cy = 



(4.55) 



4 — Differential Equations 



22 



Assume a Frobenius solutions about x = 



k+s 



and substitute into (4.55). Could you use an ordinary Taylor series instead? Yes, the point X = is 
not a singular point at all, but it is just as easy (and more systematic and less prone to error) to use 
the same method in all cases. 

oo oo oo 

(1 - x 2 ) J2 a k(k + s)(k + s- l)x k+s ~ 2 - 2x J2 a k(k + s)x k+s ~ 1 + C a kX k+s = 



oo 

S2 a k(k + s)(k + s - l)x k+s ~ 2 + 
o 

oo oo 

Y j a k [~ 2 ( k + S ) ~ (k + s){k + s- l)]x k+s + CY J < l kX k+s = 



oo 

a n+ 2{n + s + 2)(n + s + 1); 

n=-2 

oo 

«n [(n + s) 2 + (n + s)) x n+s + C^a n x 



)x n+s - 



n+s 



n=0 



n=0 



In the last equation you see the usual substitution k = n + 2 for the first sum and k = n for the rest. 
That makes the exponents match across the equation. In the process, I simplified some of the algebraic 
expressions. 

The indicial equation comes from the n = — 2 term, which appears only once. 



aos(s — 1) = 0, 



so 



0, 1 



Now set the coefficient of x n+s to zero, and solve for a n +2 in terms of a n . Also note that s is a 
non-negative integer, which says that the solution is non-singular at x = 0, consistent with the fact 
that zero is a regular point of the differential equation. 



a 2 — a 



Q"n+2 — 0"n 

s(s + l)-C 



(n + s)(n + s + 1) - C 
(n + s + 2)(n + s + 1) 

(s + 2)(s + 3) -C 



then 



O.4 — 0,2 ~ 



etc. 



(4.56) 



(4.57) 



(s + 2)(s + l)' """" ~* ~' (s + 4)(s + 3) 

This looks messier than it is. Notice that the only combination of indices that shows up is n + s. The 
index s is or 1, and n is an even number, so the combination n + s covers the non-negative integers: 
0, 1, 2, . . . 

The two solutions to the Legendre differential equation come from the two cases, s = 0, 1. 







a 



1 



-C 



x 



-C 



2-3-C 
4-3 

-c 

2 



x*+ 



2-3-C\ /4-5-C 



1 : a'n 



x + 



1-2-C 
3-2 



x 3 + 



4-3 

l-2-C\ /3-4-C 



6-5 



x 



(4.58) 



3-2 



5-4 



x 5 + ■ 
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CL2 — Qq , Q4 — (3,2 — 0, flg — Cl8 — • • • — 



and the general solution is a sum of these. 

This procedure gives both solutions to the differential equation, one with even powers and one 
with odd powers. Both are infinite series and are called Legendre Functions. An important point about 
both of them is that they blow up as x — > ±1. This fact shouldn't be too surprising, because the 
differential equation (4.55) has a singular point there. 

?t n 
y" - 7 T7 xv' + 7 V7 -,y = (4.59) 

y (l+x)(l-x) y (l + x)(l-x) y v ; 

It's a regular singular point, but it is still singular. A detailed calculation in the next section shows that 
these solutions behave as ln(l — x) near x = 1. 

There is an exception! If the constant C is for example C = 6, then with s = the equations 
(4.57) are 

-6 6-6 
2"' fl 4-a 2 — 

The infinite series terminates in a polynomial 

do + a^x 2 = ao[l — 3x 2 ] 
This (after a conventional rearrangement) is a Legendre Polynomial, 

P2(x) = -X 2 - - 

The numerator in Eq. (4.56) for a n+ 2 is [(n + s)(n + s + 1) — C], If this happen to equal zero 
for some value of n = N, then ajv+2 = ar| d so t nen all the rest of • • are zero too. The series 

is a polynomial. This will happen only for special values of C, such as the value C = 6 above. The 
values of C that have this special property are 

C = £(£ + l), for £ = 0, 1, 2, ... (4.60) 

This may be easier to see in the explicit representation, Eq. (4.58). When a numerator equals zero, 
all the rest that follow are zero too. When C = £(£ + 1) for even I, the first series terminates in a 
polynomial. Similarly for odd £ the second series is a polynomial. These are the Legendre polynomials, 
denoted Pi{x), and the conventional normalization is to require that their value at X = 1 is one. 



P (x) = 1 P 1 (x) = x P 2 (x) = \x 2 ~\ 
P z (x) = \x i -\x P^(x) = ^ x A - ^ x 2 + 



(4.61) 



The special case for which the series terminates in a polynomial is by far the most commonly used 
solution to Legendre's equation. You seldom encounter the general solutions as in Eq. (4.58). 
A few properties of the P/> are 



nm 



1 if n = m 
if n ^ m 



f 1 2 

(a) / dx P n (x)P m (x) = 2n + 1 $ 

nm where 8. 

(b) (n + l)P n +i(x) = (2n + l)xP n (x) - nP n -x{x) 

(c) Pn{x) = w^ (1 " x2)n ^ 

(d) Pn(l) = 1 Pn(-X) = (-l) n P n (x) 

00 

(e) (i-2tx + t 2 y 1/2 = Y J t n p n 



.Ax) 

n=0 
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4.12 Asymptotic Behavior 

This is a slightly technical subject, but it will come up occasionally in electromagnetism when you dig 
into the details of boundary value problems. It will come up in quantum mechanics when you solve some 
of the standard eigenvalue problems that you face near the beginning of the subject. If you haven't 
come to these yet then you can skip this part for now. 

You solve a differential equation using a Frobenius series and now you need to know something 
about the solution. In particular, how does the solution behave for large values of the argument? All 
you have in front of you is an infinite series, and it isn't obvious how it will behave far away from the 
origin. In the line just after Eq. (4.59) it says that these Legendre functions behave as ln(l — x). How 
can you tell this from the series in Eq. (4.58)? 

There is a theorem that addresses this. Take two functions described by two series: 

oo oo 
f(x) = a k xk and 9(X) = ^ b k X>t 

It does not matter where the sums start because you are concerned just with the large values of k. The 
lower limit could as easily be —14 or +27 with no change in the result. The ratio test, Eq. (2.8), will 
determine the radius of convergence of these series, and 



a k+l xk+l 



a^x k 



< C < 1 for large enough k 



is enough to insure convergence. The largest x for which this holds defines the radius of convergence, 
maybe 1, maybe oo. . . . Call it R. 

Assume that (after some value of k) all the a& and are positive, then look at the ratio of the 

ratios, 

a k+l/ a k 
bk+i/bk 

If this approaches one, that will tell you only that the radii of convergence of the two series are the 
same. If it approaches one very fast, and if either one of the functions goes to infinity as x approaches 
the radius of convergence, then it says that the asymptotic behaviors of the functions defined by the 
series are the same. 

If ? fc+1/ !? fe — 1 — ► as fast as i and if either f(x) or g(x) — > oo as x ->■ R 

°k+i/bk k 

f( x ) 

Then ——— — > a constant as x — > R 

g(x) 

There are more general ways to state this, but this handles most cases of interest. 
Compare these series near x = 1. 



— = J> & , or \ n (l-x) = J2^. 



or 



{l-x)-^ = Y j a{a - l) -- h i f t - k + l \ -x) k (a = -1/2) 



k 

k=0 
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Even in the third case, the signs of the terms are the same after a while, so this is relevant to 
the current discussion. The ratio of ratios for the first and second series is 

Qfc+i/Qfc = 1 = 1 = 1 _ 1_ 
h+i/h (k + l)/k l + l/k k 

These series behave differently as x approaches the radius of convergence (x — > 1). But you knew that. 
The point is to compare an unknown series to a known one. 

Applying this theorem requires some fussy attention to detail. You must make sure that the 
indices in one series correspond exactly to the indices in the other. Take the Legendre series, Eq. (4.56) 
and compare it to a logarithmic series. Choose s = to be specific; then only even powers of a; appear. 
That means that I want to compare it to a series with even powers, and with radius of convergence 
= 1. First try a binomial series such as for (1 — x 2 ) a , but that doesn't work. See for yourself. The 
logarithm ln(l — x 2 ) turns out to be right. From Eq. (4.56) and from the logarithm series, 

£ , s n ■, {n + s){n + s + 1) - C 

f(x)= } a n x with a n+2 = a n - — ■ — - — ■ — — 

^— ' (n + s + 2)(n + s + 1) 

n even v ' v ' 

oo 2k 

g{x) = - ln(l - x 2 ) = \ = Yl b k x2k 
To make the indices match, let n = 2k in the second series. 

E ^ = Y, Cnxn 



Now look at the ratios. 



a n+2 = n(n + l)-C = l + ^- ^ = i 2 
a n (n + 2)(n + l) 1 + ^ + ^ n 
c n+2 _ 2/(n + 2) _ n _ 1 2 



c n 2/n n + 2 l + | n 



+ 



These agree to order l/n, so the ratio of the ratios differs from one only in order l/n 2 , satisfying the 
requirements of the test. This says that the Legendre functions (the ones where the series does not 
terminate in a polynomial) are logarithmically infinite near x = ±1. It's a mild infinity, but it is still an 
infinity. Is this bad? Not by itself, after all the electric potential of a line charge has a logarithm of r in 
it. Singular solutions aren't necessarily wrong, it just means that you have to look closely at how you 
are using them. 



Exercises 

1 What algebraic equations do you have to solve to find the solutions of these differential equations? 



d 3 x dx d 10 z 
_ + a _ +te = , — 
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2 These equations are separable, as in section 4.7. Separate them and integrate, solving for the 
dependent variable, with one arbitrary constant. 

— --XN —-a 2 + x 2 — - -a(l-e-P Vx ) 
dt ~ A ' dt ~ + ' dt 1 j 



3 From Eq. (4.40) and (4.41) what are the formulas for putting capacitors or inductors in series and 
parallel? 



4 — Differential Equations 



27 



Problems 

4.1 If the equilibrium position x = for Eq. (4.4) is unstable instead of stable, this reverses the sign in 
front of k. Solve the problem that led to Eq. (4.10) under these circumstances. That is, the damping 
constant is b as before, and the initial conditions are x(0) = and v x (0) = Vq. What is the small time 
and what is the large time behavior? 

Ans: (2mv /Vb 2 + 4km)e- bt / 2m smh (^/b 2 / 'Am + k ]mt\ 

4.2 In the damped harmonic oscillator problem, Eq. (4.4), suppose that the damping term is an anti- 
damping term. It has the sign opposite to the one that I used (+bdx/dt). Solve the problem with the 
initial condition x(0) = and v x (0) = Vq and describe the resulting behavior. 

Ans: (2mvo / v '4km - b 2 )e bt / 2m sin [y/4km -b 2 tj 2m) 

4.3 A point mass m moves in one dimension under the influence of a force F x that has a potential 
energy V{x). Recall that the relation between these is 




Take the specific potential energy V(x) = —Voa 2 /(a 2 + x 2 ), where Vn is positive. Sketch V . Write 
the equation F x = ma x . There is an equilibrium point at x = 0, and if the motion is over only small 
distances you can do a power series expansion of F x about x = 0. What is the differential equation 
now? Keep just the lowest order non-vanishing term in the expansion for the force and solve that 
equation subject to the initial conditions that at time t = 0, x(0) = Xq and v x (0) = 0. 
How does the graph of V change as you vary a from small to large values and how does this same 
change in a affect the behavior of your solution? Ans: ui = \f 2Vq / ma 2 

4.4 The same as the preceding problem except that the potential energy function is +Voa 2 / (a 2 + x 2 ) . 
Ans: x(t) = Xo cosh (y / 2Vo/ma 2 1) (\x\ < a/4 or so, depending on the accuracy you want.) 

4.5 For the case of the undamped harmonic oscillator and the force Eq. (4.13), start from the beginning 
and derive the solution subject to the initial conditions that the initial position is zero and the initial 
velocity is zero. At the end, compare your result to the result of Eq. (4.15) to see if they agree where 
they should agree. 

4.6 Check the dimensions in the result for the forced oscillator, Eq. (4.15). 

4.7 Fill in the missing steps in the derivation of Eq. (4.15). 

4.8 For the undamped harmonic oscillator apply an extra oscillating force so that the equation to solve 
is 

d 2 x 

m ~dt 2 = ~ kx + F ^ t ) 

where the external force is F ext (t) = Focosuit. Assume that cu ^ coo = y/k/m. 
Find the general solution to the homogeneous part of this problem. 

Find a solution for the inhomogeneous case. You can readily guess what sort of function will give you 
a cos cut from a combination of x and its second derivative. 

Add these and apply the initial conditions that at time t = the mass is at rest at the origin. Be sure 
to check your results for plausibility: 0) dimensions; 1) oj = 0; 2) uo — > oo; 3) t small (not zero). In 
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each case explain why the result is as it should be. 
Ans: (F /m)[— cosuot + cos ut] / (cu 2 — u 2 ) 

4.9 In the preceding problem I specified that u ^ = ^Jkjm. Having solved it, you know why this 
condition is needed. Now take the final result of that problem, including the initial conditions, and take 
the limit as u — > LOq. [What is the definition of a derivative?] You did draw a graph of your result didn't 
you? Ans: (Fo/2mu)o)t sinco^ 

4.10 Show explicitly that you can write the solution Eq. (4.7) in any of several equivalent ways, 

Ae iu t + Be -iuot = Q cog UQt + D gin UQt = E cos ( Uo t + 0) 

I.e. , given A and B, what are C and D, what are E and 0? Are there any restrictions in any of these 
cases? 

4.11 In the damped harmonic oscillator, you can have the special (critical damping) case for which 
b 2 = 4km and for which u' = 0. Use a series expansion to take the limit of Eq. (4.10) as tu' — > 0. 
Also graph this solution. What would happen if you took the same limit in Eqs. (4.8) and (4.9), before 
using the initial conditions? 

4.12 (a) In the limiting solution for the forced oscillator, Eq. (4.16), what is the nature of the result 
for small time? Expand the solution through order t 2 and understand what you get. Be careful to be 
consistent in keeping terms to the same order in t. 

(b) Part (a) involved letting (3 be very large, then examining the behavior for small t. Now reverse the 
order: What is the first non-vanishing order in t that you will get if you go back to Eq. (4.13), expand 
that to first non-vanishing order in time, use that for the external force in Eq. (4.12), and find x(t) for 
small t. Recall that in this example x(0) = and x(0) = 0, so you can solve for x(0) and then for 
x(0). The two behaviors are very different. 

4.13 The undamped harmonic oscillator equation is d 2 x/dt 2 + u 2 x = 0. Solve this by Frobenius series 
expansion about t = 0. 

4.14 Check the algebra in the derivation of the n = Bessel equation. Explicitly verify that the general 
expression for a 2 k in terms of do is correct, Eq. (4.22). 

4.15 Work out the Frobenius series solution to the Bessel equation for the n = s = — l ji case. 
Graph both solutions, this one and Eq. (4.23). 

4.16 Derive the Frobenius series solution to the Bessel equation for the value of n = 1. Show that this 
method doesn't yield a second solution for this case either. 

4.17 Try using a Frobenius series method on y" + y/x 3 = around x = 0. 

4.18 Solve by Frobenius series x 2 u" + Axu' + (x 2 + 2)u = 0. You should be able to recognize the 
resulting series (after a little manipulation). 

4.19 The harmonic oscillator equation, d 2 y/dx 2 + k 2 y = 0, is easy in terms of the variable x. What 
is this equation if you change variables to z = l/x, getting an equation in such things as d 2 y/dz 2 . 
What sort of singularity does this equation have at z = 0? And of course, write down the answer for 
y(z) to see what this sort of singularity can lead to. Graph it. 
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4.20 Solve by Frobenius series solution about x = 0: y" + xy = 0. 
Ans: 1 - (x 3 /3!) + (1 ■ 4x 6 / 6! ) - (1 ■ 4 ■ 7a; 9 /9!) H is one. 

4.21 From the differential equation d 2 u/dx 2 = —u, finish the derivation for d as in Eq. (4.29). Derive 
identities for the functions c(x + y) and s(x + y). 

4.22 The chain rule lets you take the derivative of the composition of two functions. The function 
inverse to s is the function / that satisfies f[s{x)) = x. Differentiate this equation with respect to x 
and derive that / satisfies df(x)/dx = l/vl — x 2 . What is the derivative of the function inverse to 

c? 

4.23 For the differential equation u" = +u (note the sign change) use the same boundary conditions 
for two independent solutions that I used in Eq. (4.28). For this new example evaluate d and s'. Does 
c 2 + s 2 have the nice property that it did in section 4.5? What about c 2 — s 2 l What are c(x + y) and 
s(x + y)l What is the derivative of the function inverse to s? to c? 

4.24 Apply the Green's function method for the force Fq(1 — e~^*) on the harmonic oscillator without 
damping. Verify that it agrees with the previously derived result, Eq. (4.15). They should match in a 
special case. 

4.25 An undamped harmonic oscillator with natural frequency Uq is at rest for time t < 0. Starting at 
time zero there is an added force Fosinuot. Use Green's functions to find the motion for time t > 0, 
and analyze the solution for both small and large time, determining if your results make sense. Compare 
the solution to problems 4.9 and 4.11. Ans: (Fo/2ma;g) [sm(uot) — Uot cos(uot)] 

4.26 Derive the Green's function analogous to Eq. (4.32) for the case that the harmonic oscillator is 
damped. 

4.27 Radioactive processes have the property that the rate of decay of nuclei is proportional to the 
number of nuclei present. That translates into the differential equation dN/dt = —XN, where A is a 
constant depending on the nucleus. At time t = there are Nq nuclei; how many are present at time 
t later? The half-life is the time in which one-half of the nuclei decay; what is that in terms of A? 
Ans: In2/A 

4.28 (a) In the preceding problem, suppose that the result of the decay is another nucleus (the 
"daughter") that is itself radioactive with its own decay constant A2- Call the first one above Ai. Write 
the differential equation for the time-derivative of the number, iV 2 of this nucleus. You note that iV 2 
will change for two reasons, so in time dt the quantity dN% has two contributions, one is the decrease 
because of the radioactivity of the daughter, the other an increase due to the decay of the parent. 
Set up the differential equation for iV 2 and you will be able to use the result of the previous problem 
as input to this; then solve the resulting differential equation for the number of daughter nuclei as a 
function of time. Assume that you started with none, iV 2 (0) = 0. 

(b) Next, the "activity" is the total number of all types of decays per time. Compute the activity and 
graph it. For the plot, assume that Ai is substantially smaller than A2 and plot the total activity as a 
function of time. Then examine the reverse case, Ai 3> A 2 
Ans: iV Ai [(2A 2 - Ai)e~^ - A 2 e" A2t ] /(A 2 - Ai) 

4.29 The "snowplow problem" was made famous by Ralph Agnew: A snowplow starts at 12:00 Noon 
in a heavy and steady snowstorm. In the first hour it goes 2 miles; in the second hour it goes lmile. 
When did the snowstorm start? Ans: 11:23 
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4.30 Verify that the equations (4.52) really do satisfy the original differential equations. 

4.31 When you use the "dry friction" model Eq. (4.2) for the harmonic oscillator, you can solve the 
problem by dividing it into pieces. Start at time t = and position x = Xq (positive). The initial 
velocity of the mass m is zero. As the mass is pulled to the left, set up the differential equation and 
solve it up to the point at which it comes to a halt. Where is that? You can take that as a new initial 
condition and solve the problem as the mass is pulled to the right until it stops. Where is that? Then 
keep repeating the process. Instead or further repetition, examine the case for which the coefficient of 
kinetic friction is small, and determine to lowest order in the coefficient of friction what is the change in 
the amplitude of oscillation up to the first stop. From that, predict what the amplitude will be after the 
mass has swung back to the original side and come to its second stop. In this small /Lt|< approximation, 
how many oscillations will it undergo until all motion stops. Let b = fMkFjy Ans: Let t n = im/uiQ, 
then for t n <t < £ n +i. x(t) = [xo — (2n + l)b/k] coswot + {—l) n b/k. Stops when t « irkxo/2uob 
roughly. 

4.32 A mass m is in an undamped one-dimensional harmonic oscillator and is at rest. A constant 
external force Fq is applied for the time interval T and is then turned off. What is the motion of the 
oscillator as a function of time for all t > 0? For what value of T is the amplitude of the motion a 
maximum after the force is turned off? For what values is the motion a minimum? Of course you need 
an explanation of why you should have been able to anticipate these two results. 

4.33 Starting from the solution Eq. (4.52) assume the initial conditions that both masses start from the 
equilibrium position and that the first mass is given an initial velocity v x \ = Urj- Find the subsequent 
motion of the two masses and analyze it. 

4.34 If there is viscous damping on the middle spring of Eqs. (4.45) so that each mass feels an extra 
force depending on their relative velocity, then these equations will be 

m 1 d ^ = -k 1 x 1 -k 3 (x 1 -x 2 )-b(x 1 -x 2 ), and 

m2 ~dW = - k 3 (x 2 - xi) - b(x 2 - xi) 

Solve these subject to the conditions that all initial velocities are zero and that the first mass is pushed 
to coordinate Xq and released. Use the same assumption as before that m\ = m 2 = m and k\ = k 2 . 

4.35 For the damped harmonic oscillator apply an extra oscillating force so that the equation to solve 
is 

dP x dx 

= -b Tt - kx + F ext (t) 

where the external force is F ext (t) = Foe tujt . 

(a) Find the general solution to the homogeneous part of this problem. 

(b) Find a solution for the inhomogeneous case. You can readily guess what sort of function will give 
you an e %Lut from a combination of X and its first two derivatives. 

This problem is easier to solve than the one using coswt, and at the end, to get the solution for the 
cosine case, all you have to do is to take the real part of your result. 



4.36 You can solve the circuit equation Eq. (4.37) more than one way. Solve it by the methods used 
earlier in this chapter. 



4 — Differential Equations 



31 



4.37 For a second order differential equation you can pick the position and the velocity any way that 
you want, and the equation then determines the acceleration. Differentiate the equation and you find 
that the third derivative is determined too. 

d 2 x b dx k . |. d 3 x b d 2 x k dx 
dt 2 m dt m dt 3 m dt 2 m dt 

Assume the initial position is zero, x(0) = and the initial velocity is ^(0) = v$\ determine the second 
derivative at time zero. Now determine the third derivative at time zero. Now differentiate the above 
equation again and determine the fourth derivative at time zero. 

From this, write down the first five terms of the power series expansion of x(t) about t = 0. 
Compare this result to the power series expansion of Eq. (4.10) to this order. 

4.38 Use the techniques of section 4.6, start from the equation md 2 x/dt 2 = F x (t) with no spring 
force or damping, (a) Find the Green's function for this problem, that is, what is the response of the 
mass to a small kick over a small time interval (the analog of Eq. (4.32))? Develop the analog of 
Eq. (4.34) for this case. Apply your result to the special case that F x (t) = Fq, a constant for time 
t > 0. 

(b) You know that the solution of this differential equation involves two integrals of F x (t) with respect 
to time, so how can this single integral do the same thing? Differentiate this Green's function integral 
(for arbitrary F x ) twice with respect to time to verify that it really gives what it's supposed to. This is 
a special case of some general results, problems 15.19 and 15.20. 

4.39 A point mass m moves in one dimension under the influence of a force F x that has a potential 
energy V(x). Recall that the relation between these is F x = —dV/dx, and take the specific potential 

111 

energy V(x) = —Vq€~ x ' a , where Vo is positive. Sketch V . Write the equation F x = ma x . There 
is an equilibrium point at X = 0, and if the motion is over only small distances you can do a power 
series expansion of F x about x = 0. What is the differential equation now? Keep just the lowest 
order non-vanishing term in the expansion for the force and solve that equation subject to the initial 
conditions that at time t = 0, x(0) = Xq and v x (0) = 0. As usual, analyze large and small a. 

4.40 Solve by Frobenius series methods 

d 2 y 2dy 1 _ 

dx 2 x dx x^ 

Ans: Er=o(- 1 ) fc ^TT)I is one - 

4.41 Find a series solution about x = for y" + ysecx = 0, at least to a few terms. 
Ans: a [l - \x 2 + Ox 4 + ^x 6 + •••]+ ai[x - |x 3 - ^x 5 + • • • ] 

4.42 Fill in the missing steps in the equations (4.55) to Eq. (4.58). 

4.43 Verify the orthogonality relation Eq. (4.62)(a) for Legendre polynomials of order £ = 0, 1, 2, 3. 

— 1/2 

4.44 Start with the function (l — 2xt + t 2 ) . Use the binomial expansion and collect terms to get 
a power series in t. The coefficients in this series are functions of x. Carry this out at least to the 
coefficient oft 3 and show that the coefficients are Legendre polynomials. This is called the generating 
function for the P e 's. It is ^ Pl{x)t l 
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4.45 In the equation of problem 4.17, make the change of independent variable x = l/z. Without 
actually carrying out the solution of the resulting equation, what can you say about solving it? 

4.46 Show that Eq. (4. 62)(c) has the correct value P n (l) = 1 for all n. Note: (l-x 2 ) = (l+x)(l-x) 
and you are examining the point x = 1. 

4.47 Solve for the complete solution of Eq. (4.55) for the case C = 0. For this, don't use series 
methods, but get the closed form solution. Ans: Atanh -1 x + B 

4.48 Derive the condition in Eq. (4.60). Which values of s correspond to which values of £? 

4.49 Start with the equation y" + P{x)y' + Q(x)y = and assume that you have found one solution: 
y = f{x). Perhaps you used series methods to find it. (a) Make the substitution y(x) = f(x)g(x) 
and deduce a differential equation for g. Let G = g' and solve the resulting first order equation for G. 
Finally write g as an integral. This is one method (not necessarily the best) to find the second solution 
to a differential equation. 

(b) Apply this result to the I = solution of Legendre's equation to find another way to solve prob- 
lem 4.47. Ans: y = f J dx -h exp — j P dx 

4.50 Treat the damped harmonic oscillator as a two-point boundary value problem. 

mi + bx + kx = 0, with x(0) = and x(T) = d 
[For this problem, if you want to set b = k = T = d = 1 that's o.k.] 

(a) Assume that m is very small. To a first approximation neglect it and solve the problem. 

(b) Since you failed to do part (a) — it blew up in your face — solve it exactly instead and examine the 
solution for very small m. Why couldn't you make the approximation of neglecting ml Draw graphs. 
Welcome to the world of boundary layer theory and singular perturbations. Ans: x(t) « e 1 ~ t — e 1 ~ t / m 

4.51 Solve the differential equation x = Ax 2 (I + tut) in closed form and compare the series expansion 
of the result to Eq. (4.25). Ans: x = a/[l- Aa(t + ut 2 /2)] 

4.52 Solve the same differential equation x = Ax 2 (l + ut) with x(to) = a by doing a few iterations 
of Eq. (4.27). 

4.53 Analyze the steady-state part of the solution Eq. (4.42). For the input potential Voe lujt , find the 
real part of the current explicitly, writing the final form as I max cos(cot — <fi). Plot I max and (j) versus 
to. Plot V(t) and I(t) on a second graph with time as the axis. Recall these V and / are the real part 
understood. 

4.54 If you have a resistor, a capacitor, and an inductor in series with an oscillating voltage source, 
what is the steady-state solution for the current? Write the final form as I max cos(ut — (p), and plot 
/ max and (j) versus cu. See what happens if you vary some of the parameters. 

Ans: / = V cos(ut-(j))/\Z\ where \Z\ = ^ R 2 + {uoL - l/uoC) 2 and tan0= (uL - l/uC)/R 

4.55 In the preceding problem, what if the voltage source is a combination of DC and AC, so it is 
V(t) = Vq + V\ cosut. Find the steady state solution now. 
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What is the total impedance left to right in the circuit ^yy^ R 2 



Ans: R 1 + (l/iuC 2 ) + l/[(l/iwLi) + l/((l/iwC , 1 ) + l/((l/i? 2 ) + (l/iwL 2 )))] 



4.57 Find a Frobenius series solution about x = for y" + ycsca; = 0, at least to a few terms. 

A nc . rtf 1 rr» 2 I 1 /y*3 1 /y*4 I 1 /y»5 I 

iio. 2 12 48 192 

4.58 Find a series solution about x = for a; 2 ?/' — 2irq/ + (x 2 + i — l)y = 0. 

4.59 Just as you can derive the properties of the circular and hyperbolic trigonometric functions from 
the differential equations that they satisfy, you can do the same for the exponential. Take the equation 
u' = u and consider that solution satisfying the boundary condition u(0) = 1. 

(a) Prove that u satisfies the identity u(x + y) = u{x)u(y). 

(b) Prove that the function inverse to u has the derivative l/x. 

4.60 Find the asymptotic behavior of the Legendre series for the s = 1 case. 

4.61 Find Frobenius series solutions for xy" + y = 0. 



